Geometric mesoscopic correlations in quasi-one dimension 
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We study analytically and numerically field/intensity correlations in wave transport through 
volume-disordered waveguide. The obtained channel and spacial correlations deviate from those 
found in framework of Dorokhov-Mello-Pereyra-Kumar (DMPK) formalism, that we relate to inap- 
plicability of equivalent channel approximation in DMPK. We show that this can be remedied by 
introducing boundary correction - an escape function which depends on the waveguide geometry - 
that describes wave transport near a boundary between random medium and free space. We obtain 
the expressions for field/intensity channel and spacial correlation functions which agree with the nu- 
merics and are consistent with the perturbative expressions in slab geometry as well as experiments 
conducted in QID. 
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Mesoscopic fluctuations, such as universal conductance 
fluctuations 1] , are rooted in nonlocal correlations 0, Q 
that appear due to the interference effects when a wave 
undergoes multiple scattering events in random medium. 
Although common to both the electron and the electro- 
magnetic waves, the experiments with the light give more 
detailed information about transport 0, II] , for example 
mesoscopic correlation of optical fields leads to readily 
observable speckles [1]. 

The universality of the statistics of wave transport and 
ability to study localization-delocalization transition as 
a function of only system length, made QID geometry 
{i.e. a wire) a fruitful testbed for studying mesoscopic 
phenomena 0, B Q ■ Due to quantization of the trans- 
verse momentum, the wave transport can be described in 
terms of Tab = ^abf^b ~ the channel transmission coeffi- 
cient - that measures the flux transmitted in the channel 
b when a unit flux enters via channel a. Summation over 
the channel indexes gives the transmission Ta ~ Y^bTab 
and the dimensionless conductance [IQ] : 



T = '^Tab^'^kakbTab; g = (T) . 



(1) 



where ka,b are the longitudinal component of the mo- 



mentum, necessary for the proper llj treatment of in- 
coming/outgoing fluxes; Tab is the squared amplitude of 
the mode b. For the light, Tab, Ta and T are directly 
measurable 13, 1^ . 

Diagramma tic |l| and Dorokhov-Mello-Pereyra-Kumar 

(DMPK) [3, [Hi techniques were highly successful in 

studying mesoscopic phenomena and were applied to in- 
vestigate the correlations in the slab [l3| and QID [l3| 
geometries respectively. DMPK is a self-embedding ap- 
proach that assumes the incremental transfer matrices to 
be random with only time reversal and energy conserva- 
tion being the constraints. Such random matrix theory 
(RMT) relies on isotropy assumption, which leads to a 
mathematically convenient anzatz [3] but lacks a micro- 
scopical foundation 18. [l^. This approximation yields 



results which depend on number of channels N, irrespec- 
tive of geometry of the waveguide. It can be recast into 
the equivalent channel approximation (EGA), validity of 
which was questioned before [H, [3, [lO] and recently 
received a renewed focus [IJl. EGA assumes that after 
scattering in a macroscopically thin slice of the disordered 
medium all channels become completely mixed and, thus, 
are statistically equivalent. This approximation is linked 
to DMPK's inability 



21\ to adequately describe dif- 
fusion in the direction transverse to the axis of the wire, 
as well as a discrepancy between mean freepath (MFP) 
i values in RMT and the transport theory [9|. |l9l|. Mod- 
ified DMPK formalism, free of EGA, yields [21| complex 
equations which resist an analytical treatment. 

The purpose of this work is the detailed test of 
RMT-KECA predictions in QID [i,[l3| with the emphasis 
on mesoscopic correlations. Using the direct numerical 
simulations, we obtain the transmission matrix without 
making EGA or any other approximations. We show that 

1/2 1/2 

in the N ^ oo limit, it is the quantity tabPa Pb that 
exhibits the properties of RMT transfer matrix: (i) its 
elements are identically distributed for different channel 
indexes a, 6; (ii) field and intensity correlation functions 
match those derived in RMT. For L ^ £, the escape 
function pa depends only on the geometry of the waveg- 
uide, independent of its length L. This suggests that 
the limitations of EGA can be remedied by the multi- 



plicative correction pa^"^ pl^"^ 



to t 



(RMT) 



We show that 



the resultant spacial correlations adequately describe the 
diffusion in the crosscction of the waveguide and are con- 
sistent with diffusion theory. 

Numerical method: In our numerical simulation, we 
consider 2D waveguide filled with random medium [i^, 
[ist . We use finite difference time domain method to cal- 
culate the response of our system to pulsed excitation, 
followed by Fourier transformation which gives us the 
desired continuous- wave response [i^l so that tab can be 
computed 1^^. Because the studied system is QID, the 
samples with values of g in the range {0.4 — 4} were ob- 



FIG. 1: (a) Mean value of the transmission normalized by 
conductance N "^^^iTab) / "^^fjiTab) obtained from numerical 
simulations in samples with different A'^, g, and £: (i) O - 
N ^ 19, g ^ 3.5; (ii) □ - iV = 19, g ^ 1.5; (iii) A - 
N = 19, g = 0.5; (iv) <) - N ^ 15, g = 0.8; (v) < - 
N = 28, g ^ 3.0; (vi) > - N = 19, g = 4.5, £ which is 
twice as long as in the previous samples. The data agrees 
with p{/Lt) and A — 0.818) (thick solid and dashed lines) 
without any adjustable parameters. Thin solid line (a con- 
stant) depicts the Equivalent Channel Approximation made 
in DMPK approach. The inset shows spacial field correla- 
tion function (see Eq. |8]) obtained in numerical simulation 
(squares), Eq. (|8} (thick line) and an expectation based on 
EGA (thin line). 



tained by varying the lengths L of the random medium. 
To ensure that the obtained results do not depend on 
the microscopic structure of disorder, we operate in the 
regime of locally weak disorder, k£ 1. To analyze 
statistics of mesoscopic transport, ensembles of 10'* — 10^ 
random realizations were obtained. 

Validity of EC A : We directly verify the earUer reports 
dl) Hj EO) El that (Tab) has strong channel dependence. 
Fig. [T] shows that it is described by the Chandrasekhar 
function as found for a slabfBI. 



{Tab) 

p (p) — C exp 



pika/k) X p{h/k) X {T)/N^ (2) 
fi ln(l -cos/?(/3/sin/3)^-2 

TT 7o cos^ /3 + sin^ (3 



where £) = 2, 3 is the dimensionahty of the waveguide, 
and C is chosen so that pa is normalized Pa = N. 
Throughout this work (...) imply averaging only over dis- 
order configurations for given channel indexes {a, b}. The 
agreement between Eq. ^ and the numerical data in 
Fig. [T]is achieved with no adjustable parameters for sam- 
ples with different N, g, and ^ with the only condition 
L ^ £ to avoid the ballistic regime The deviations 
from p in Fig. [1] are attributed to limited statistics and 
finite number of channels N . 

Such significant channel dispersion should not come 
as a surprise as it reflects fundamental wave coherence 
properties [IJl , and has a long history in radiative trans- 



FIG. 2: Statistical equivalence of the normalized transmission 
coefficients. Depicted are two groups of a = [l..Af] histograms 
P{sab = Tab/ {Tab}) (b is sct to 1 for example) and P{sa = 
Ta/{Ta)) for the sample with g = 3.5, = 19. Both groups 
are successfully fitted by the RMT's expression from Ref. |24l ]. 



fer theory (RTT) [28|. For slab geometry p{fi) oc p + A 
is known to be a good approximation with A — zt/i is 
the so-called extrapolation length. A equals to 0.818 and 
0.710 (in 2D and 3D respectively) arise in RTT, whereas 
Ad = 7r/4 « 0.785 and 2/3 w 0.667 are obtained in the 
diffusion approximation [5|. Also, the persistent channel 
dependence of (Ta) was seen in the simulations of Ref. 
[29|, but was interpreted as a finite N effect. 

Statistics of transmission coefficients: RMT's relation 
{Tab) = {T)/N'^ stems from the assumption of the sta- 
tistical equivalence of all channels - EGA. In view of Eq. 

this connection is no longer a trivial one. E.g. in 
the surface disordered waveguides, strong channel de- 
pendence of MFP can even lead to coexistence of dif- 
fusive and localized regimes [30| in different channels 
of the same sample. In contrast, in case of volume- 
disorder, considered here, the disorder-induced channel 
mixing should lead to equilibration for L ^ £. Inde- 
pendence of Pa from g and L [l9| may be considered a 
manifestation of this effect. 

Our simulations. Fig. [2l demonstrate that all channels 
indeed become statistically equivalent when the trans- 
mission coefficients are normalized by their averages - 
P{sab = Tab/{Tab)) and P{sa = Ta/{Ta)) arc indepen- 
dent of channel indexes. It suggests that removing pa.b 
factor is sufficient to enforce EGA. It is, thus, tempting 
to associate tab/ Pa"^ / p^"^ or Sab with t^ab^^'^^ from RMT. 
However, Sab are no longer required to obey the compo- 
sition rules 



(3) 



The discrepancy appears due to correlations in different 
channels, to be discussed below. The knowledge of Eq. 
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A + 2/tt ^ 
A + 7r/4 " 
A + 1/2 ^ 
A + 2/3 " 


i 0.907 
i 0.879 


Ad + 2/tt ^ 
Ad + 7r/4 " 
Ad + 1/2 _^ 
Ad + 2/3 " 


a 0.905 
i 0.875 



TABLE I: a and /3 defined in Eqs. (|4ll0p quantify tlie devia- 
tion from EGA. Tlie values in tfiree columns were computed 
in A'^ ^ oo limit, where '-^'^ replaced with the appro- 
priate integration. p{fi) from Eq. @ was used in the first 
column. Approximation p(/i, A) with A from RTT and dif- 
fusion theory were used to compute the values in the second 
and the third columns respectively. 
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FIG. 3; Ghannel intensity correlations i!TabTa'b') normalized 
by ^ah) and {Ta'b') to remove dependence on p^.t. General 
structure of the resultant correlator remains the same in all 
studied samples, regardless of g. Sample with g = 0.4, N — 
19 is shown in the figure. To visualize the four-dimensional 
array a' = 7, b' = 12 was chosen. 



^ is sufficient to derive the following scaling relations 




1 



N 
(4) 



where a = (l/N) is a measure of non-equivalence 

of the channels. Using pa or pa gives a > 1, Table HI As 
one can see, the breakdown of the composition rules Eq. 
([3]) is in itself a mesoscopic effect related to the channel 
dispersion. Eq. ^ also shows that the normalization in 
the definitions of Sab and Sa is sufficient to remove the 
channel dependence in studies of mesoscopic fluctuations 
in ^ cxD limit. 



Channel- channel correlations: Diagrammatic pertur- 
bative calculations of the channel correlations in the slab 






FIG. 4: Two types of channel correlators and their leading 
diagrams in l/g perturbation series that appear in derivation 
of spacial intensity correlation function, Eq. ([9]). The known 
expression for (TafcTk'b') is shown in the first line, see e.g. 
Leading contributions to the correlator {tabtab''^a'b''ta'b) '^an 
be obtained by swapping the incoming and outgoing indexes 
as shown in the second row, Eq. (O follows. 

geometry d, [^, [l^ give a general relation {TabTa'b') « 
{Tab) {Ta'b') ■ Here a,b denote different transverse mo- 
menta. The averages {Tab) appear when the incoming 
and outgoing paths are paired into diffusons Q , the only 
contributions that survive averaging over disorder. The 
fact that {Tab) has nontrivial dependence on a, h in QID 
has not been widely appreciated. It is the lack of this de- 



pendence that makes (nonperturbative) RMT result [17 1 
not directly applicable. Comparison with the results of 
the previous section shows that modifications needed to 
correct this shortcoming may not be straightforward. 

We find channel-resolved field and the intensity corre- 
lations in the form (Fig. [3|) 

{tabt*aib') — ~ — '^^j2^^ 5aa'5bb'\ {TabTa'b') — (5) 



{T)'^ PaPbPa' Pb' 

7V4 



iV2 

X [ASaa'Sbb' + B{Saa' + Sbb') + A] (6) 



where A,B are some functions of g. The structure of 



the above expression coincides with RMT prediction [17 1 
only in EGA, i.e. by setting pa = 1. Therefore, multi- 
plication by PaPbPa'Pb' is Sufficient to recover the proper 
correlations function, which does describe the diffusion 
within the crossection of the waveguide. It is also con- 
sistent with the expressions obtained with diagrammatic 
techniques in the slab geometry [1, [^, [l^ where pa,b be- 
come continuous function of p. Upon multiplication by 
Pa,bS, the correlations Eqs. (|5l6p regain the dependence 
on the waveguide dimensionality and its shape (in 3D) 
that is in contrast to RMT where only number of chan- 
nels enters as a sole parameter. 

Spacial correlations: Spacial correlations of the trans- 
mitted fields and intensities in QID have been a focus of 
the experimental and numerical studies [13, [US [3]. 
The results were compared to diagrammatic calcula- 
tions in slabs. Despite the agreement with the the- 
ory, this may not be well justified because in QID 
(where the waveguide width W is smaller then L) the 
boundary conditions in the transverse direction(s) dif- 
fer from those in a slab. At the same time, starting 
with the RMT's (which enforces proper boundary condi- 
tions through quantization of the transverse momentum) 



4 



expectation {tabtl,i,,) = {{T) /N'^)Saa'Sbb' , accounting for 
the variation of the propagation speed for different modes 
and then averaging over the output cross-section of the 
waveguide gives the foUowing expression (2D) Ce{^v) — 
{Eiyo + Ay)E*{yo)) cx i:^^,il/h) cosiirbAy/W). It 
does not agree neither with our simulation nor diffusion 
theory [3l|, inset in Fig. [U and leads to unphysical non- 
zero correlations oc 1/N for large separations. 

We will see that account of the channel depen- 
dence in Pa's corrects the discrepancy. In case 
of intensity correlations, the starting point in mak- 
ing transformation from channel to spacial coor- 
dinates is the correlator {taibitl^b^jasbatlibj- % 
pairing the channel indexes into diffuson contri- 



butions one can see that this correlator contains 

four terms 5aia2^bib2^a3a4^f)3&4 + ^aia4^f)i&4^a2a3^62b3 + 
'5aia2^blb4<^a3a4^b2b3 + <^aia4 <5bib2 '^0203 <Jb3b4 • The firSt tWO 

lead to Eq. ([6]), whereas the last two require the knowl- 
edge of the new correlator (iab^ab'^a'b'^a'b), which accord- 
ing to our diagrammatic analysis, Fig. Ul has form: 

{{Tf/N^) PaPbPa'Pb' X [BSaa'Sbb' + A{Saa' + 6bb') + B] 

(7) 

with the same A, B as in Eq. ([6]). This is fully supported 
by our simulations for all studied samples. Based on 
Eqs. (|5I6I7P we derive the following expressions for the 
spacial field and intensity correlation functions in a 2D 
waveguide: 



{EiYa + Ay, j/o + Ay)E* (Fq, yo)) = 



(/(ro + Ay,yo + Ay)/(yo,2/o)) = 



f{AY) /(Ay); /(y) = 1 ^ ^ cos a y 



X 



4 {A fiAY) f{Ay) + B [f{AY) + f{Ay)\ + A] 
I 



(8) 
(9) 



where 7[pa] = TliaPo-l^^o- ^'2/ ^'^'^ respectively the 
transverse position of the source and detector in the per- 
fect lead regions before and after the disordered part of 
the waveguide. In the above, averaging over Iqi 2/o where 
performed and only 0(1/A^") terms were retained. The 
derived Eq. ([5]) gives an excellent fit to data, inset in 
Fig. [T] Furthermore, substituting p(/x, A^ and changing 
summation over p^ to the integral over p-JA gives ex- 
actly the diffusion expression for 2D slabs [23|, |3l| . The 
structure of Eq. ([9]) also agrees with the perturbative 
expression [s^ where A = 1 + 2/15g^ and B = IjZg. 

Setting Ay and Ay to zero in Eqs. (|8l9p gives the 
first two moments of the transmitted intensity. Their 
analysis shows that channel dependence in pa also affects 
the actual value of the average intensity. Indeed, RMT's 
expression not corrected by Pa would be reduced by a 
factor 



l\Pa - 1] 



\ ^ Pa_ I \ ' 1 



(10) 



which is computed Table U as /3 < 1. 

Conclusion: The following two observations provide 
an insight into why simple multiplicative correction by 
Pa was able to successfully remedy the limitation of 
RMT set by EGA: (i) The function pa is independent 
of system size, and thus is insensitive of the nature 
of wave transport inside random medium; (ii) spacial 
field correlation function predicted by RMT in 2D (the 
same conclusion holds also in 3D) gives Bessel function 
CsiAy) cx S^i(l/fcb)cos(7r6Ay/IF) = Jo(fcy) which co- 
incides with Shapiro [35| expression for the correlations 



inside random medium. Therefore, we conclude that pa^b 
describes the surface effect^ whereas RMT properly de- 
scribes the transport in the bulk of the random medium. 
Multiplication of RMT's transfer matrix by pa,b ~ the 
escape function - makes the geometry dependent bound- 
ary correction that leads to correct result in iV — s- cxd 
limit. We derived Eqs. (|5I6I7I8I9P based on the RMT 
with the surface correction that is consistent with known 
expressions for slab geometry. The obtained expressions 
also agree with our direct numerical simulations in both 
y < 1 and y > 1 regimes. Albeit the studied correla- 
tions are directly measurable only in optics, our results 
are also applicable to noninteracting electronic systems - 
a volume-disordered wire, where time reversal symmetry 
is preserved and dephasing mechanisms can be neglected. 
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